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ABSTRACT 



For single-file translocations of a charged macromolecule through a narrow pore, the cru- 
cial step of arrival of an end at the pore suffers from free energy barriers, arising from changes 
in intrachain electrostatic interaction, distribution of ionic clouds and solvent molecules, and 
conformational entropy of the chain. All contributing factors to the barrier in the initial stage 
of translocation are evaluated by using the self-consistent field theory for the polyelectrolyte 
and the coupled Poisson-Boltzmann description for ions, without radial symmetry. The bar- 
rier is found to be essentially entropic, due to conformational changes. For moderate and 
high salt concentrations, the barriers for the polyelectrolyte chain are quantitatively equiv- 
alent to that of uncharged self-avoiding walks. Electrostatic effects are shown to increase 
the free energy barriers, but only slightly. The degree of ionization, electrostatic interaction 
strength, decreasing salt concentration and the solvent quality all result in increases in the 
barrier. 
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I. INTRODUCTION 

of single polyelectrolyte molecules 
through narrow pores is one of the most fundamental processes encountered in many bi- 
ological processesi"^ and technological^'^'^'^'^'^'^»^»^»^ applications. An inevitable need to 
understand the translocation phenomenon has driven the scientific community to study 
both natural^i'^i^'^i^'^i^ as well as synthetio^'i2i»ii'i^ polyelectrolytes. In these studies, an ex- 
ternal driving force is used to carry out a successful polyelectrolyte translocation, which 
is typically due to an applied electric Heldr-^^^^^^^^'^^^^ and in some cases, arise due to 
an osmotic imbalance (confinement-driven translocatioi>i'2i»^»^»^»^»^). In addition to the 
complications coming from an intricate coupling between the short range excluded volume 
interactions and the long-range electrostatics, a satisfactory description of the translocation 
of a polyelectrolyte chain must take into account many extra factors. In general, translo- 
cation of a polyelectrolyte chain from one confined space to the other through a narrow 
pore may be affected by the dielectric mismatch between the interior and exterior of the 
confining membrane^»iiii^i2i^ nature of the por^'^^^^^, surface of the confining membrane^, 
electro-osmotic fiow through the transmembrane pore^^, and semi-fiexibility of the chain^^. 
Furthermore, depending on the area of cross-section of the pore, the polyelectrolyte may 
undergo translocation either as a single-file with linear conformations or as multiply folded 
conformation. 

Independent of the nature of the driving force and any additional factor mentioned above 
that may affect translocation, the single-file translocation is envisioned as a two step process-. 
In the first step, one end of the chain arrives at the pore entrance^i^i^i^ and, in the second 
step, the chain is th.Tea.ded:^^^^^'^^ through the pore from one side to the other. It has 
been recognized that both steps are associated with entropic barriers, with the first barrier 
associated with the loss of translational entropy of chain ends and the second barrier with 
reduction in conformational entropy of the chain. 

The extent of the entropic barrier due to conformational changes among various con- 
tributing factors for the experimentally relevant polyelectrolytes is not known. In fact, it is 
known^^ that the conformational entropy of a confined polyelectrolyte is only a weak contrib- 
utor to the free energy where translational entropy of small ions and molecules (counterions, 
colons, and solvent) is dominant. Reorganization of counterion clouds around deforming 
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polyelectrolyte chains can also contribute significantly to the free energy. Furthermore, the 
intrachain electrostatic repulsion can stiffen the polymer enabling an easier access to the 
pore entrance. It is therefore of interest to assess the relative magnitudes of various con- 
tributing factors to the free energy barrier associated with the translocation of a flexible 
polyelectrolyte molecule. 

Although most of the theoretical works on the entropic barrier for translocation has fo- 
cused on the threading (second step), the largest part of the barrier is actually associated 
with the first step of localizing one end of the chain at the pore entrance^.. In this work, 
we focus on the first step in confinement- driven translocation involving single flexible poly- 
electrolyte chain trying to get out of the confining spherical cavity through a pore on the 
surface as a "single-file" and provide a quantitative description of the free energy barriers 
for the chain end to find the pore. The translocation barrier for the first step is estimated 
by the free energy difference between "end-fixed" (one end fixed near the pore) and "free 
ends" (confined chain, which is free to move inside the cavity) equilibrium states of the chain 
(Fig. [1]). Using this approach, the free energy barriers for a Gaussian chaii>i^»^»^, trapped 
initially inside a spherical cavity, can be computed exactly and are purely entropic in nature 
due to lower degrees of conformational freedom in "end-fixed" state as compared with "free 
ends" and the absence of interactions. Similar calculations for the excluded volume chain^^ 
have been carried out within spherical symmetry and in the absence of solvent. For the 
case of polyelectrolyte translocation, these barriers are unknown and form the focus of this 
study. 

Earlier theories of translocation of polyelectrolytes have been constructed only by using 
results of neutral polymer chainsF^'^^ and ignoring the coupling between conformations of 
the polyelectrolyte chain and the small ions. However, it is widely being recognized that the 
physics of polyelectrolytes is dominated by counterions. In view of this, it becomes necessary 
to estimate the role of counterions and small electrolyte ions in establishing the free energy 
barriers for translocation. Here, we present a systematic calculation of free energy barrier by 
an explicit treatment of the coupling between small ions and conformations of polyelectrolyte 
chains. By adopting the self-consistent-field theory (SCFT) for a fiexible polyelectrolyte 
chain and combining with the Poisson-Boltzmann prescription for the electrolyte ions and 
counterions, we have computed the various energetic and entropic contributions to the free 
energy barrier. Since the localization of one chain end at a specific location on the surface of 
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the cavity breaks the radial symmetry, we have solved the self-consistent coupled nonlinear 
differential equations in two dimensions with azimuthal symmetry. As pointed out earlier, we 
address only the localization of one of the chain ends at the pore, without any consideration 
of all effects arising from the pore itself. 

The rest of the paper is organized as follows: theory is presented in Sec. Ull numerical 
technique is presented in Sec. Illlj calculated results and conclusions are presented in Sec. 
|IV]and lYl respectively. 

II. THEORY 

In order to study confinement driven translocation, we consider a single negatively 
charged flexible polyelectrolyte chain in a spherical cavity of radius R and model the chain 
as a continuous curve of length Nb, where is the number of Kuhn segments, each of length 
b. An arc length variable t is used to represent any segment along the chain backbone so that 
t G [0, A^]. To maintain global electroneutrality, we assume that the spherical cavity is filled 
with ric monovalent counterions (positively charged) released by the chain in addition to 
ions of species 7 (= +, — ) coming from added salt. Moreover, we assume that there are 
solvent molecules (satisfying the incompressibility constraint after assuming the small ions 
to be pointlike) present in the cavity of volume Q and for simplicity, each solvent molecule 
occupies a volume (vs) same as that of the monomer (i.e., Vs = b^)- Subscripts p, s, c, + and 
— are used to represent monomers, solvent molecules, counterions from the polyelectrolyte, 
positive and negative salt ions, respectively. The valency (with sign) of the charged species 
of type j is represented by Zj and the degree of ionization of the chain is taken to be a. 
In the point charge limit for the small ions considered here, the cations from the added 
salt (j = + ) and the counterions of the polyelectrolyte chain (j = c) are indistinguishable 
from each other. Also, we consider smeared charge distribution so that each of the segments 
carries a charge eaZp, where e is the electronic charge. 

We use self-consistent field theory (SCFT) to compute the free energy of single flexible 
polyelectrolyte^^ chain in "free ends" and "end-fixed" states. Ignoring the potential inter- 
actions between solvent molecules and small ions, the partition function for the single chain 
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system in either of the states can be written as 
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where R(t) represents the position vector for t*^ segment and subscripts = s,c, +,— . 
What distinguishes the "end-fixed" state from the "free ends" state in Eq. ( [1]) is the 
functional integral over R. Physically, the functional integral over R represents the sum 
over all the possible conformations of the chain originating from one end and ending at the 
other. Explicitly, for the "free ends" state, the functional integral is given by J -D[R] = 
J dvQ J drj^ J^^ -D[R], where fq and r^v are the positions of the ends of the chain represented 
by the specific values of the contour variable t = and t = N along the chain, respectively. 
Similarly, /D[R] = J dr^ J^^ D[Ii] for the chain, whose one end is fixed at fq in the 
"end-fixed" state. 

Note that in Eq. ( [I]), it is understood that the factor of 1/2 in the last term inside 
the exponent is present, only when j = j' and ksT is the Boltzmann constant times abso- 
lute temperature. Furthermore, V^p(r), ^^^(r) and l^s(r) represent the interaction energies 
for monomer- monomer, solvent-solvent and monomer-solvent pairs, respectively, when the 
interacting species are separated by distance r =| r | and are given by 

P ^ (2) 



Vpp( 


» 


= Wpp6{r) + 






= WssSir), 


Vpsi 




= Wps5{r). 



AneoekBT r 

(3) 
(4) 

In writing the interaction energies, the short range excluded volume interactions are mod- 
eled by three dimensional delta functions [5(r)] multiplied by the respective excluded volume 
parameters. For the monomer-monomer, solvent-solvent and monomer-solvent pairs, these 
parameters are taken to be Wpp, Wss and Wps, respectively. Also, the long range electrostatic 
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interactions are modeled by Coulomb's law after assuming the effective dielectric constant 
(e) of the medium to be position independent, eo being the permittivity of the vacuum. 

In writing Eq. ( [1]), the constraint on the number densities of the monomers and the 
solvent molecules to obey the incompressibility condition at all points inside the spherical 
cavity is written as the product of delta functions involving microscopic densities on the right 
hand side (po being the total number density of the system so that po = {N + ns)/Q = 
) . As mentioned earlier, the incompressibility condition is written after taking the small ions 
to be point-like. For the point-like limit of the small ions, the interaction energies between 
the monomers and the ions, represented by Vpj, are given by 

Vp,{r) = for j = c,+,-. 5 

Similarly, the interaction energies between the small ions can be written as 

^"'W-JS^f 

Following the protocol presented in Appendix A, the free energy F of the single chain 
system in either of the states can be computed using the well-known saddle point approx- 
imation. F is expressed as integrals over inhomogeneous number densities of the various 
components of the system and the electric potential. Taking the dielectric constant (e) of 
the medium to be independent of temperature (T) and Flory's chi parameter defined as 
Xpsb^ = — {wpp + Wss) /'^ ~ 1/7", the free energy (within saddle-point approximation) can 
be divided2^(see Appendix B and C for details) into enthalpic part due to excluded volume 
and electrostatic interactions, and entropic part due to small ions, solvent molecules and 
the polyelectrolyte chain. Denoting these contributions by E^, E^, S ions , S solvent and Spoiyi 
respectively, the free energy can be written as 

Eu] -\- Eg T (^Sions ~l~ S solvent ~l~ Spoly)i (T) 

where Fq = ^ {Nwpp + UsWss)- Superscript -k denotes that the saddle point approximation 
has been used to compute the free energy. Explicit expressions for the different contributions 
are given by 

Evj = Xpsb^ j drpp{r)ps{r), (8) 
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\ J dr^(r)pe(r), (9) 

-r5.o;.ent = y" drp,{r) {In [p,(r)] - 1} , (11) 

-TSpoiy ^ -]nQp- Jdr [{^^^^(r) + Wp{r)} pp{r)] . (12) 

In these equations, Pd{r) and Wisi^r) are respectively the macroscopic number density and 
the field experienced by species of type f3, due to excluded volume interactions at the saddle 
point. Also, all the charged species experience an electrostatic potential represented by ip{r), 
which is related to the local charge density Pe(r) = J2j=c + - ZjPj{r)+Zpapp{r) by Poisson's 
equation, 

V^V'(r) = -^TThPeir). (13) 

Note that ■ip{r) in these equations is dimensionless (in units of kBT/e) and Ib is the 
Bjerrum length defined as = /4TreoekBT. 

At the saddle point, the macroscopic densities for the small molecules are related to the 
corresponding fields by the Boltzmann law, so that 

= ".exp|-u.(r)] 
J drex.p[-Ws{r)] 

J drexp[-ZjV(r)] 

The fields and densities are related to each other by the saddle point equations, given by 

Wp{r) = Xpsb^Ps{r)+viT^), (16) 
Ws{r) = Xpsb'^Ppir) +r]{r), (17) 

where 77(r) is the Lagrange's multiplier introduced to enforce the incompressibility constraint. 
For the "free ends" state, the monomer density is dependent on the field by the relation 



drq{r, N) 

and for the "end-fixed" state, the relation becomes 

-TV 



(18) 
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Superscript / and a depict the free and anchored nature of the single chain. These monomer 
densities are related to the solvent density by the incompressibility constraint 

Pp{r) + ps{r) = po. (20) 

In Eqs. ( US]) and ( [12]), the function q{r,t) is the probability of finding segment t at 
location r, when starting end of the chain can be anywhere inside the spherical cavity, and 
it satisfies the modified diffusion equation 



dq{r,t) 



,2 



b' 



g(r,t), (21) 



dt 

along with the initial condition g(r, 0) = 1. Similar to q{r,t), the Green function, 
G(r, ra,t, 0), is the probability of finding segment t at location r, when starting end of the 
chain is at ra. It also satisfies Eq. ([21]) but with the initial condition G(r, Fa, 0, 0) = 5(r— ra), 
where S represents the three dimensional delta function. The partition function of the chain 
(Qp) can be written in terms of these functions. Specifically, for the "free ends" state, it 
is given by Qp = = J drq{r, N) and for the "end-fixed" state, the partition function 
becomes Qp = = J drG{r, r^, N, 0). 

We include the effect of confinement by solving Eqs. ([I3]- [2T|) using the Dirichlet bound- 
ary conditions for Ws{r),ip{r), q{r,t) and ^(r. Fa, t, 0). Rationale for choosing these bound- 
ary conditions is that we model the spherical cavity as a neutral, hard surface. Dirichlet 
boundary conditions for tp{r) corresponds to a neutral spherical cavity. Similarly, Dirichlet 
boundary conditions for q{r,t) and G(r, Tg, t,0) means the spherical cavity is like a hard 
wall and the monomer density at the surface must be zero. This is equivalent to carrying 
out calculations with a wall potential, which is a delta function. The reason for choosing 
the Dirichlet boundary condition for Ws{r) is the following: physically, the incompressibility 
constraint gets violated near the surface of the confining spherical cavity and due to the 
hard wall like spherical cavity, monomer as well solvent density at the surface must be zero. 
Masking techniques'^ have been used to take care of this fine point and study the surface 
effects. In this work, we take a different approach. Noticing that the free energy contribu- 
tion from the values of densities and fields at the hard surface is zero, we choose Ws{R) = 
and compute the solvent density at the surface. Note that pinning Ws{R) to zero also means 
that the field is shifted by a constant. However, it is well-known that the densities and the 
free energy at the saddle-point are independent of this shift. 
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For the "free ends" state, these equations have been solved after assuming radial 
symmetry^. Since the radial symmetry is broken when one end of the chain is anchored 
near the surface at the site of the pore, we have solved the SCFT equations in two dimen- 
sions with the azimuthal symmetry. In two dimensions, due to inherent singularities at the 
center of the spherical cavity (due to the division by zero) in addition to those arising from 
the delta functional form for the initial condition of the modified diffusion equation for the 
"end-fixed" state, the numerical solution of SCFT equations is non-trivial. In the present 
work, we have developed an efficient numerical scheme to solve these coupled non-linear 
equations in the presence of aforementioned singularities, which is presented in the next 
section. 



III. NUMERICAL TECHNIQUE 

As mentioned earlier, we need to solve the above non-linear set of equations in two dimen- 
sions under the assumption of azimuthal symmetry. Note that the "end-fixed" state has an 
azimuthal symmetry in a polar coordinate system, where the origin is at the intersection of 
the surface of the confining sphere and the radius connecting the center of the sphere to the 
anchoring point. This means that in the coordinate system, where the origin is the center of 
the confining sphere (Fig. Wp), the SCFT equations need to be solved on a semi-circle with 
its diameter along the x-axis. The solution on a semi-circle can be obtained by defining 6 as 
the angle with respect to the only diameter of the semi-circle so that the SCFT equations 
need to be solved in {r, 9} space, where r G [0,R],9 E [0, n]. In practice, this procedure can 
be carried out using = ^ [^('^^^) + ii^M i^^^^^)]' "^^^^ means that in order to use 
the pseudo-spectral method^li^ for an accurate solution of the modified diffusion equation, 
we need to use the Legendre polynomials and hence, the Legendre transform. This partic- 
ular step slows down the computational procedure due to the unavailability of an efficient 
fast Legendre transform. 

In order to develop a faster algorithm, we note that for the particular value of 6 = 7r/2, 



the Laplacian becomes V 



2 _ J_ 
r 



Due to the particular functional form 



of the (j) dependent part in the Laplacian, fast Fourier transform^' (FFT) can be used 
to apply this part of the Laplacian. In other words, we need to solve on a circle so that 
r G [0, i?], G [0, 2ti\ rather than a semi-circle. However, the use of FFT in the computations 
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over the Legendre transform speeds up the calculations even when the number of collocation 
points on the grid gets enlarged by a factor of four. Also, note that we expect the results 
(for the densities etc. ) on a circle to be symmetric about the axis passing through the 
anchoring point and the center of the sphere due to the azimuthal symmetry of the problem. 
This provides a nice check on the computational procedure. 

In view of this computational efficiency, we have solved the SCFT equations using spher- 
ical polar co-ordinates {r,6 = 7r/2,0) so that r G [0, -R] and (p G [0,27r]. Instead of solving 
these equations in real or Fourier space, we use the split-step pseudo-spectral method^^ii^^ 
employing the fast Fourier transforn>^ (FFT) and sine transforn>2£, which allows a faster 
and accurate computation of the densities and free energies. It is convenient to solve for 
f{r,t) = rq{r,t) rather than solving for q{r,t) directly. As earlier, the solution of Eq. ( [2T|) 



is obtained after using 



dr 



dr ' 



and expanding /(r,t) = Y.i=Qfi{"^^'^y ■> 



where L is the number of terms required to represent the function / by a finite series within 
the desired accuracy. Writing Eq. ( [2T|) in terms of f{r,t) and using Baker-Hausdorff 
formula^^i^^, the solution is given by the propagation relation 



f{r,t + dt) 



exp [—dtw{r)/2] exp 
62 92 



exp 



dt- 



6 dr'^ 



exp 



dt 
dt 



12 r2 
62 1 92 



12 r2 



exp [-dtw{r)/2]f{v,t), 



(22) 



where w{y) = Zpaipijc) + Wp{v). Also, the multiphcation by r leads to f{0,(f),t) = for all 
values of t so that numerical problems at r = 0, due to division by in the Laplacian are 
avoided. Another major advantage of the transformation is that now, the equations are to 
be solved with periodic boundary conditions. So, FFT and sine transform can be used to 
implement the exponential of operators. Exponential of (p dependent operator on the right 
hand side of Eq. ( [22|l is applied in Fourier space after taking one forward and one backward 
FFT. r dependent operator is applied after taking forward sine transform defined as 

K 

fi{r,t) = J29k{l,t)sm{knr/R), (23) 



k=l 



and the corresponding inverse sine transform, so that boundary conditions for / are always 
satisfied during the computations. Here, K is the number of terms in the finite series 
on the right hand side of Eq. ( [23]) to accurately sample the function fi{r,t) within the 
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desired accuracy. In total, each time step requires two FFTs with respect to and one 
sine transform with respect to r in each direction (forward and backward). The same 
technique has been used for solving G{r, r^, t, 0) after approximating the initial condition 
for ^(r, Fa, 0, 0) = (5 (r - r^) by 

5(r-ra) = I ^^rf^ (24) 
[O if r^r,. 

Here, A(j) = 2tt/L and Ar = R/K are the grid spacings used to discretize the two dimensional 
space. 

To solve Poisson's equation, Vr'?/'(r) = — 47r/BPe(i"), we use a similar strategy. We solve 
for h{r) = rip{r), so that Poisson's equation becomes 

hir) = -AnlBrpeir). (25) 
Now, expanding h{r) = Xljlo ^K'")^*''^' ^^e equation for components hi becomes 

hi{r) = -47r/BFFT<^ [rpe(r)] . (26) 

Here, the subscript means, FFT is to be taken with respect to 0. These sets of equations 
can be readily solved for the real and imaginary parts of hi{r) due to the tridiagonal nature 
of the finite difference equation set obtained with the constraints hi{0) = hi{R) = 0. Now, 
taking backward FFT with respect to 0, h{r) is obtained. 

Starting from an initial guess for fields, Wp{r),Ws{r) and '0(r), new fields and densities 
are computed using the method described above with the boundary conditions mentioned 
in the previous section. Simple mixing-^ is used to obtain the new guess and the iter- 
ative process is continued until the difference in newly computed and the guessed fields 
is of the order of 10~^. Using the converged solution for the fields and densities, free 
energies for the "free ends" (i.e., F* = F^) and "end-fixed" {F* = F*) states are com- 
puted. For the computation of the free energy barriers for the chain end to find the 
pore, we compute the difference AF = F* — F^. To analyze different contributions to the 
free energy barriers, we also compute the differences in the enthalpic and entropic contri- 
butions, given by AE^ = ( "end-fixed" ) - E^("free ends"), A^e = Se( "end-fixed") - 
Ee( "free ends"), -TAS'io„^ = -rS'io„^( "end-fixed" )+ rS'io„^( "free ends"), -TAS'pojj^ = 
-7'5'poZjy( "end-fixed" )+ TS'poiy( "free ends") and -TASsoivent = -T SsoiventC' end- fixed" ) + 
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T5'so«?;ent("free ends"). The results presented here were obtained by using L = 32, K = 32 
and dt = 0.1 after optimizing the numerical algorithm for speed and accuracy. Also, we have 
fixed one end of the chain for "end-fixed" state at ra = [(i? — 0.625)6, 7r/2, 0] in spherical 
polar co-ordinates . 

IV. RESULTS 

In Sec. ini we have presented the theoretical analysis for a flexible polyelectrolyte chain 
in the presence of counterions and colons of arbitray valency. However, it is well-known in 
the literature that the saddle-point approximation used in this work (which provides the 
Poisson-Boltzmann description for the electrolytes) fails^i^i^^ in the case of multivalent 
ions. For monovalent ions'^^, Poisson-Boltzmann equation provides quite reasonable results. 
In view of this, we have considered here a polyelectrolyte chain with monovalent counterions 
in the presence of monovalent salt. 

For monovalent counterions and colons, we have carried out extensive numerical compu- 
tations of the free energy barriers by varying the various parameters required to solve the 
above coupled equations namely, A^, R, a, Ib, Xps, and the salt concentration Cg (in units 
of moles per liter, i.e., M). For wide ranges of a, Ib, Xps, and Cg, we have solved the above 
equations for the A^— and R— dependencies of the monomer, counterion and colon density 
profiles, electric potential distribution and the various free energy contributions. It turns 
out that these extensive numerical calculations lead to some general conclusions, presented 
below. 

A. Monomer, counterion and coion distributions 

Typical monomer and electrostatic potential distributions for the "free ends" and "end- 
fixed" states are shown in Fig. [2] for A^ = 50 and R/b = 4. Here, we have chosen 
a = 0.1,1b = ^b,Xps = 0.45, and Cg = O.IM, where the choice of the parameters has been 
motivated by the experimental relevance to aqueous systems. Although these parameters 
cannot be varied independently in experimental situations, we have computed the conse- 
quences of each of these parameters in order to obtain physical insight into the origins of 
the free energy barriers. It is evident from Figs. [2^ and that, for the unanchored state, 



13 



that both distributions are radially symmetric, with the monomer density and the negative 
electric potential being a maximum at the center of the cavity. As one end of the chain is 
localized near the right-edge of the equator (represented by the arrow in Figs. [Jt and [2]i), 
the density and potential distributions become anisotropic as expected. 

Also, the coupling between the counterion, colon and monomer distributions, whose origin 
lies in the electrostatic interactions between these charged species, can be seen clearly in 
Fig. [31 In the figure, we have plotted the counterion and colon distributions for the same 
single chain systems, whose monomer and potential distributions are shown in Fig. [21 For 
the "free ends" state, it is found that the counterion and colon distributions (Figs. [3^ 
and [3)d, respectively) are also radially symmetric like the monomer distribution shown in 
Fig. [3^. Furthermore, Figs. [2^ and [3^ show that the counterion distribution tracks 
the monomer density distribution with the maximum number density of the counterions at 
the center of the cavity. On the other hand. Fig. [3)d shows that the colons are excluded 
from the regions rich in monomer density with a minimum number density of the colons at 
the center of the cavity. These results are consistent with the calculations carried out for 
the "free ends" state with radial symmetry^. However, Figs. [3t and [3li show that the 
radial symmetry gets broken for the "end-fixed" state and the anisotropy of the small ion 
distributions follow from the electrostatic coupling between counterion, colon and monomer 
distributions. In addition, it is to be noted that although the counterions and colons are 
distributed differently for the "free ends" and "end-fixed" states, the net electric potential 
( Figs. [2}d and [2li) tracks the monomer density distribution. 

These results (Figs. [2l and [3]) show that the anchoring of the chain end at a specified 
location leads to anisotropic monomer and charge distributions. Intuitively, we anticipate 
that larger anisotropics in the monomer and charge distributions for the "end-fixed" state 
relative to the "free ends" state correspond to larger free energy differences between the states 
and hence, larger free energy barriers. Also, due to the coupling between the monomer and 
small ion distributions, a longer chain in "end-fixed" state is expected to show a higher 
degree of anisotropy in the monomer and charge distributions compared to a shorter chain. 
However, in a confined space, the increase in the chain length leads to space filling and 
counteracts the anisotropic distributions, which is otherwise expected. We demonstrate this 
particular point in Fig. [Hby presenting the monomer and electrostatic potential distribution 
for the chain with twice the number of monomers as in Fig. [2l (i.e., = 100) and keeping 
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all other parameters the same. 

Comparing Figs. [2^ with H^, and ^ with Ht, it is clear that the increase in the 
chain length for a fixed radius of the confining spherical cavity leads to an increase in 
the monomer density at all points in the interior for both the "free ends" and "end-fixed" 
states as expected. Furthermore, monomer density distributions for "free ends" and "end- 
fixed" states (Figs. H^i and Ht, respectively) are found to be almost indistinguishable from 
each other except a small region near the anchoring point. Similar trends are seen for the 
electrostatic potential (compare Fig. HJd with Sli) and small ion distributions (not shown). 

To be more quantitative, we have plotted the monomer density profiles in the "free ends" 
and "end-fixed" states along x and y axes for different values of for a given spherical 
cavity (Fig. [5]). We have plotted these density profiles for the values of at which we see 
noticeable differences in the free energy barriers for the polyelectrolyte and the self-avoiding 
walk chains. From Figs. and \5jp, it is clear that the system has radial symmetry in "free 
ends" state. However, the radial symmetry is broken in "end-fixed" state as it is clear from 
Figs. [St and [5li. Also, as mentioned earlier the monomer density increases everywhere 
inside the spherical cavity, when A^ is increased for both the "free-ends" and "end-fixed" 
states. 

Furthermore, the expected increase in the degree of anisotropy for the monomer and 
charge distributions in the "end-fixed" state with the increase in the chain length does not 
occur for all monomer volume fractions in a confined system, where space filling effects coun- 
teract the anisotropy. This particular point is demonstrated in Fig. [6], where we have plotted 
the monomer density distributions along x and y axes in the "free ends" and "end-fixed" 
states at different monomer volume fractions. At a low volume fraction (corresponding to 
A^ = 50 at = 4 in Figs. Efe- and [Hb), there is significant anisotropy in the monomer 
density distributions arising from the fixing of one end of the chain. However, at a higher 
volume fraction (corresponding to A^ = 100 at = 4 in Figs. and Eli), the monomer 
density distributions in the "free ends" and "end-fixed" states become almost indistinguish- 
able from each other. Note that the electrostatic potential and the density distribution of 
the small ions also become indistinguishable from each other in the two states at higher 
volume fractions (Figs. [7^ and [Tb, respectively). This is a result of the space filling or 
confinement effect counteracting the anisotropic effects originating due to the localization of 
one of the ends of the chain. This particular point will be used later to explain the trends 
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in the free energy differences in the "free ends" and "end-fixed" states presented below. 

B. Free energy barriers 

In Fig. [HI we have plotted the free energy difference (in units of ksT) between the "end- 
fixed" and "free ends" states for different values of and R. For comparison purposes, 
the barriers for the chains when electrostatics is switched off (i.e., self-avoiding walk chain 
with a = 0, Cs = 0, Ib = 0), are also plotted. It is evident that the free energy barriers 
for polyelectrolyte chains are almost identical to those for the corresponding uncharged self- 
avoiding walk chain at higher monomer densities and they differ only by a small amount in 
lower density regime. As seen in Fig. [HI the dependence of the free energy barrier on the 
chain length is nonmonotonic, for a given cavity size. 

The nonmonotonicity in the free energy barriers has also been seen in the case of self- 
avoiding walk chains in the absence of sol vent In agreement with Ref.— , the origin of 
the nonmonotonicity lies in the entropic, excluded volume and the confinement effects. For 
low enough values of such that the net intra-chain excluded volume interaction is weak, 
the probability^ of finding a particular monomer (such as the chain end) at any prescribed 
spatial location decreases with an increase in A^. This is equivalent to an increase in the free 
energy barrier with A^ in this limit. The entropic contributions to the free energy barriers 
arising from the anisotropic small ion distribution add to this effect (see the description 
below on the effect of electrostatics in Figs. [TOl[TT]and [T2|) . On the other hand, for higher 
packing fractions of monomers, the excluded volume and the synergistic space filling effects 
coming from the confinement take over and consequently the two equilibrium states become 
indistinguishable from each other (compare Figs. [2] and [4]). This is equivalent to a decrease 
in the free energy barrier with an increase in A^ for the limit of strong confinement. 

Although the numerical values of and F* are significantly different for the polyelec- 
trolyte and uncharged polymer cases, the agreement between the free energy barriers for 
these two cases, as seen in Fig. [HI is striking. In order to identify the origins of the barriers 
for polyelectrolytes and of the agreement with uncharged polymers, we have plotted different 
energetic and entropic contributions for two radii in Fig. [91 It is found that the dominant 
contribution to the barriers is the difference in conformational entropy {—TASpoiy) of the 
chain in the "free ends" and "end-fixed" states. Difference in solvent entropy {—T /S.S solvent) 
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and energy (AE^) due to excluded volume interactions between different constituents plays 
a role, although meager, only when chain length is small and solvent is the major component 
in the system. Contributions due to the difference in entropy {—TASions) of small ions and 
electrostatic energy (AEe) are negligible in comparison with other contributions. 

It has already been shown that the excluded volumer interaction energy (Ey^), electrostatic 
energy {Eg) and conformational entropy {—TSpoiy) of a flexible polyelectrolyte chain are 
minor contributors to the absolute chain free energy^. Only for very strong confinements 
(with monomer volume fractions higher than 0.7), the chain conformational entropy starts 
contributing significantly. Otherwise, it is the entropy of the small ions ( — TSions) snd 
the solvent molecules {—TS solvent), which dominate the free energy. We have confirmed 
the same results using the two dimensional algorithm presented in Sec. IIIII (comparison 
not shown here). Based on the relative contributions to the total free energy coming from 
different components in the system, one might think that the change in counterion, colon or 
solvent distribution would dominate the free energy barriers. However, on the contrary, we 
have found that the dominant contribution to the free energy barriers is the change in the 
conformational entropy of the chain. 

Analysis of the different contributions to the free energies, calculated for the differ- 
ent values of the parameters, leads to the following explanation. Entropies of small ions 
{—TSions) and the solvent molecules {—TSgoivent) are strongly dependent on the number of 
these molecules in the system and show a very weak dependence on the spatial distributions 
(cf. Eqs. ( [TOj) and ( [TT]) ). For a single chain system in the presence of salt, a large number 
of small ions and solvent molecules explain the dominance of the total free energy by their 
entropies. However, due to an equal number of small ions and solvent molecules in the "free 
ends" and "end-fixed" states, the entropies of the small ions and the solvent molecules are 
almost the same in either of the states and cancel each other almost exactly in the compu- 
tation of the free energy barriers. This explains the minor contribution to the free energy 
barriers due to the small ions and solvent molecules in Fig. [91 

Similarly, from Eq. ( [8]), the excluded volume interaction energy (Eui) depends on the 
monomer distribution. However, Figs. [2^, and Ht reveal that an increase in the 

chain length for a given radius of the cavity makes the "end-fixed" state almost indistin- 
guishable from the "free ends" state in terms of monomer distribution due to confinement 
effects. Consequently the excluded volume interaction energy contributions to the free en- 
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ergy barriers are very small at higher monomer volume fractions. Numerical results also 
show that the electrostatic energy (Ee) contribution to the total free energy is orders of 
magnitude lower than the other contributions (e.g., Ee is of the order of 0.1 for = —400 
for R/b = 4 in Fig. [8]). Keeping in mind such a low contribution to the total free energy, 
the negligible contribution to the free energy barriers (AF) from the change in electrostatic 
energy (AEf.) is not a surprise. However, the chain conformational entropy depends on the 
distribution of the chain ends (cf. Eq. [T2]) and indeed differ in the two states due to a 
relatively lower number of conformational states available for the chain in the "end-fixed" 
state compared to the "free ends" state. As a result, the difference {—TASpoiy) shows up 
as the dominant contribution to the free energy barriers (AF) in Fig. [91 Furthermore, 
comparing Figs. [9^ and [9)d, it is found that in addition to the dependence of the barriers 
directly on the monomer density, the cavity radius plays an important role in modifying the 
barriers. In agreement with the work on Gaussian chaina^^, the cavity radius provides an 
entropic contribution to the barriers going like 41n(i?/Z). This explains the decrease in the 
free energy barriers with a decrease in R, as is evident in Fig. [81 

These results indeed support the entropic nature of the free energy barriers. Furthermore, 
these results are quite robust (within 1 — 2A;bT) for a vast range of parameters involving 
a, Is/b, Xps and c^. In Fig. [TUl we have plotted the free energy barriers for different values of 
the degree of ionization (a) of the polyelectrolyte chain. It is found that indeed electrostatics 
play a role in the free energy barriers and that the barriers increase with the increase in the 
degree of ionization of the chain at lower volume fractions. However, the increase in the free 
energy barriers is small (e.g., the change in the free energy barriers is less than 10% when 
a is changed from 0.1 to 0.5 in Fig. [TOl) . Reason for the small change in the absolute value 
of the free energy barriers is the already mentioned logarithmic dependence of the small ion 
entropy term (i.e., —TASions) on the anisotropic ion distributions (cf. Eq. [TD|) and a weak 
dependence of the chain entropy {i.e.,—TASpoiy) on a. Although the absolute change in the 
free energy barriers is small when the degree of ionization is changed by a factor of 5, it is 
worthwhile to investigate the origin of the change in the free energy barriers. It is found that 
the change in the free energy barriers has a dominant (but small in magnitude) contribution 
coming from the small ion entropy (e.g., —TASions changes from 0.03 to 0.44 compared to 
the change of —TASpoiy from 6.11 to 6.28 when a is changed from 0.1 to 0.5 in Fig. [TUlat the 
lowest volume fraction). In other words, the slight increase in the free energy barriers comes 
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mainly from the anisotropic distribution of the small ions in the "end-fixed" state. At higher 
monomer volume fractions, the confinement effects take over and the barriers are the same 
as for the self-avoiding walks. Similarly, the effect of added salt concentration and Bjerrum 
length on the free energy barriers are shown in Fig. [TT] and [121 respectively. Similar to the 
effect of a, the free energy barriers increase with the lowering of the salt concentration and 
an increase in Bjerrum length i.e., with the strengthening of the electrostatics. Note that 
the change in the free energy barriers due to the change in Bjerrum length is minuscule due 
to the weak contribution of the electrostatic energy (cf. Eq. [HD to the free energies^i. 

The effect of solvent quality on the free energy barriers is shown in Fig. [131 From the 
figure, it is clear that an increase in the solvent quality i.e., a decrease in Xps, leads to an 
increase in the free energy barriers. This is an outcome of the change in excluded volume 
interaction energy with the change in the solvent quality (see Eq. [HD and can be explained 
as following. With the change in the solvent quality, it is found that all the contributions 
to the free energy barriers remain almost the same except the excluded volume interaction 
energy (i.e., AE^ given by Eq. [HI). Note that the excluded volume interaction energy 
depends quadratically on the monomer density distribution (because Xpsb^ J 'j^rpp(r)ps(r) = 
Xpsb^N — Xpsb^ J drpp{r) as a result of the incompressibility constraint). Due to the higher 
monomer density near the surface of the confining cavity in the "end-fixed" state compared 
to the "free ends" state, AE^ is negative (cf. Fig. [HD- Also, the prefactor Xps in Eq. [HI 
causes a decrease in the magnitude of A^^^^, (which is negative) with the decrease in Xps (i-e., 
the increase in the solvent quality). In view of this, the free energy barriers increase with 
the increase in the solvent quality. 



V. CONCLUSIONS 



In summary, we find the remarkable result that even though the conformational entropy of 
a flexible polyelectrolyte chain is a minor contributor to the chain free energy, the free energy 
barrier is essentially entirely due to the change in the conformational entropy of the chain, 
for experimentally relevant conditions of translocation experiments. Even more remarkably, 
the free energy barrier for a flexible polyelectrolyte for moderate salt concentrations is not 
significantly different from that for an uncharged self-avoiding chain. The free energy barrier 
increases with degree of ionization, Bjerrum length, and solvent quality, and decreases with 
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salt concentration. However, the increase in the free energy barriers in the confined single 
chain system investigated here is small. Furthermore, it is to be noted that the entropically 
driven free energy barrier for placing one chain end at the pore entrance is about 6 — 9 ksT, 
which is within the access of energy released in one event of ATP hydrolysis^. Nevertheless, 
in experiments involving fast translocations, the barriers computed here with equilibrium 
consideration might be modified by non-equilibrium polymer conformations. 

Finally, it must be remarked that the present development of the model and numerical 
scheme to treat the anisotropy of polyelectrolyte conformations is only a starting point. The 
influences of electrostatic forces arising from dielectric mismatches due to the pore-bearing 
membrane, and the nature of the charged pore itself are some of the future directions for 
integrating theories of polyelectrolyte translocation into experimental investigations. 
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APPENDIX A : Self Consistent Field Theory 

Here, we present the details about the calculation of the free energy of a single poly- 
electrolyte chain within saddle point approximation. Similar procedure has been used 
earlier— i^i^i^'^. In order to carry out the transformation from a description involv- 
ing particles to the fields, we define a dimensionless Flory's chi parameter as Xpsb^ = 
u!ps — {wpp + Wss)/2 along with microscopic densities as 

Pp{v) = / dt5[v-n{t)], (A-i) 

Jo 

pA^) = ^^("^ -^i) fo^^ j = -5,C, +,-, (A-2) 

i=l 

Pe(r) = aZpPp(r) + ^ ^iPj(r)> (A-3) 

i=c,+,- 

where pp{r), pj{r) and epe(r) stand for monomer, small molecules (ions and solvent 
molecules) and local charge density, respectively (in units of e, e being the charge of an 
electron) . 
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Using these definitions, the partition function in Eq. ( [T]) can be rewritten as 



exp 



F~Fo 



D\K] / n n 

j m=l 

I 



exp <; / dt 



N 



'dTj 



-Xpsb J drpp{r)ps{r) - y / dr / dr 
Y[6 {pp{r) + ps{r) -po) , 



,Pe(r)pe(r') 



(A-4) 



where Fo/ksT = (po/2) {Nwpp + risiUss) and Ib = /AireQekBT. 

To carry out the transformation from particles to fields, we use the following three trans- 
formations in order: (1) the well-known Hubbard-Stratonovich transformation for the elec- 
trostatics part, which leads to the introduction of field ip in the calculations by 
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(2) functional integral representation for unity to decouple the excluded volume interac- 
tions between the monomers and solvent molecules 
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k=p,s 



n ^Nexp 

k=PyS 



J d^Y^ Wk{r) {pfc(r) - pfc(r)} 

k=p,s 



(A-7) 



which leads to the introduction of collective fields Wk and densities pk for the monomer and 
solvent molecules. Also, the transformation leads to the replacement of microscopic density 
variables (p^) by the collective density variables (pfc). 

(3) functional integral representation for the delta functions to enforce incompressibility 
constraint at all points in the system 



l[5 {pp{r) + p^iv) 



Po) 



D[ri{r)] exp 



drvir) (pp(r) +p,(r) - po) 



(A-8) 



which introduces the well-known pressure field ri{r) in the calculations. 

Using these transformations along with Stirling's approximation In n! ~ n In n — n, Eq. 

( IA-4p becomes 
p p ^ 

^ h — rt c k=p,S 
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U D[wk] j n D[p,] I D[r]] j D[tP]exp[-f{wk,Pk,V,^}], 
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where 



f {wk,pk,V,i'} = Xpsb^ j drpp{r)ps{r) +i j dr7]{r) {pp{r) + p,{r) - po) - InQp 
-i dr^ Wk{r)pk{r) + ^ 



k=p,s 



j=s,c,+,- 



^3 _1 



In 



1 



Stt/ 



(ir-(/'(r)Vr'?/'(r), 



(A-10) 



where Qp, Qs and Qj are the single particle paritition functions for the polyelectrolyte chain, 
solvent molecule and small ions of different species. Explicitly, the chain partition function 
is given by 

"N To / \ 2 



L'[R(t)]exp 



(A-11) 



Similarly, single particle partition function for a solvent and small ion is given by 

Qs = j dr exp[-iws{r)], (A-12) 



J dr exp [—iZjip(r)] for j = c, +, 



(A-13) 



The functional integrations over the collective fields and densities in Eq. ( IA-9P are almost 
impossible to compute exactly. A well-known approximation to evaluate the functional 
integrals is the steepest descent technique (also known as saddle-point approximation). We 
use the approximation to compute the free energy of the single chain system so that the 
approximated free energy is given by F — ~ — Fq = / {w^, pl,r]*,il!*} (after taking 
ksT = 1 as a scale for energy), where collective fields and densities with stars as superscripts 
are their respective values at the saddle point and lead to the extremum of the functional 
/. Carrying out extremization of the functional / with respect to the collective fields and 
densities, Eqs. ( [131- EQ]) are obtained and details of carrying out the standard functional 
derivativesS^iSI. Note that the normalization factor p^ is ignored in the calculations for the 
approximate free energy by the saddle-point approximation. 



APPENDIX B : Energy and entropy of single polyelectrolyte chain within saddle 

point approximation 

Here, we present the calculations of the energetic and entropic contributions to the free 
energy of a single polyelectrolytic chain using thermodynamic arguments. The analysis is 
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similar to the one presented in Ref.— , where polymer conformations and hence, the monomer 
density are kept fixed during the well-known Debye charging process. For the free energy 
within the saddle point approximation used in this work, polymer conformations and hence, 
the monomer density are also dependent on the charging parameter and need to be treated 
properly. 

Computation of the energetic and entropic contributions to the final free energy within the 
saddle point approximation can be carried out after imagining two isothermal "charging" 
processes. One is the traditional electric charging process, where charges of all the ions 
in the system (free ions in the solution as well as on the chain backbone) are increased 
gradually from to their final values. Similarly, we imagine another "charging" process, 
where excluded volume parameter characterizing the short range excluded volume effects is 
developed (or "charged") gradually from to its final value, which, in turn, leads to the 
development of monomer density and field waves. 

For the electric charging, we assume that at any instance, charge of the ions is a fraction 
A of its final value. Similarly, the excluded volume "charging" process is characterized by the 
charging parameter A, so that at any instance the excluded volume parameter characterizing 
the excluded volume interactions between species i and j is w[j = X^Wij, where Wij is the 
final value of the excluded volume parameter. Here, we follow the notation used by Marcus^^ 
to represent quantities at any instance in the charging process by superscript ' and the choice 
of A^ to relate the instant value of the excluded volume parameter to the final value is made 
to simplify the mathematics as discussed at the appropriate places in this Appendix. Also, 
the order of charging for the two processes does not matter as expected. 

Before computing the energy and entropy for the single polyelectrolyte chain with the 
solvent treated by using the incompressibility condition, it is worth investigating the chain 
in the absence of solvent so that there is only one excluded volume parameter Wpp charac- 
terizing the monomer-monomer excluded volume interactions. In this particular case, the 
field experienced by monomers is Wp{r) = WppPpir). Note that here Wpp has the dimen- 
sions of volume^^ and is inversely proportional to the temperature. If the local charge per 
unit volume at any instance during the electric charging process is represented by p'^^iy) = 
A y^,- ^ , Zip'Ar)\ ^ + (r)| , , , where 4j'(r) and w'(r) = w'p'(r)\ , , , 

are the electrostatic potential and the field arising from excluded volume interactions at the 
particular instance, then work done [dFe) in charging a small volume element dQ = dr by 
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amount dX is given by dF^ = dXdr [Ei=c,+,- ZjPji^)\x^'ir) + ^p^Ppi^)\x^'(r),^'^ir)\ ksT^p'ir) = 
^drp'^{r)kBTtlj'{r). So, total work done in charging the whole system from A = to A = 1 
is given by 

where 

drp'e{r)i^'{r) (B-2) 

is the electrostatic energy at any instance during the electric charging process. 

Similarly, the work done in development of excluded volume parameter d\^ for a small 
volume element dfl is given by dFyj — dr^Pp{r))\^^,^^^ ^, ^^^kBTWp{r) , ks and T being Boltz- 
mann's constant and temperature, respectively. Hence, the total work done in developing 
the final value of Wpp for the whole system (i.e., from A = OtoA = l)is given by 



F^ = / d\'^ (B-3) 
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keT 

is the energetic contribution to the free energy at any particular instance arising from the 
excluded volume interactions. 

From thermodyamics, we know that the free energy {F) at any temperature T is related 
to energy {E) and entropy (S) by the following relations: 

F = E-TS (B-5) 
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Integrating the last relation between the energy and the free energy, 
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F = T / Ed[-] (B-7) 
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It can be shown^^ that the work done as calculated from the isothermal charging processes 
(i.e., Fe + F^) is related to the free energy (F) calculated by integrating energy E over the 
inverse temperature (cf. Eq. IB-7|1 and in general, Fg + Fy^ is not equal to F. For the two 



processes considered here, the total work done during isothermal charging (Few) is given by 
Few = Fe + Fw Few cau bc related to F by noting that the expression for the energy at any 
instance is the same in the two descriptions i.e., E' = E'^ + E'^, because from Eq. ( IB-7p . 



d{F'/T) 
dl/T 

Also, from Eqs. ([B3]) and f lB^ . 
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dX^ J J. X^T 
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Here, F^ and F^ are the work done during the electric and excluded volume charging at any 
instant characterized by a particular value of the charging parameter A. 
From Eqs. i^^, f lB^ and f lB^ . 



2 f diF'/T) \ f d{Fi + Fj2)/T \ 



It can be shown that Eq. ( IB-121) is satisfied as long as (Fg + F^/2)/T is related 
to A and T only through A^/T at fixed dielectric constant e (after noting that ^ = 
J^^^'^^ d ^In ^ j ^). Marcus has already shown that the use of the Poisson-Boltzmann 
equation for computing the electrostatic potential does not lead to the violation of the con- 
straint presented in Eq. ( IB-12p . Also, this particular constraint for relating the work 
done from the isothermal charging process to the free energy has led us to choose A^ as the 
prefactor for instant value of the excluded volume parameter so that the instant values of 
the monomer densities and fields are dependent on A^/T (after taking the excluded volume 
parameter Wpp to be of the form a/T, a being a constant independent of temperature T). 
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Furthermore, note that if the instant value of excluded volume parameter is taken to be 
= Xwpp similar to the electric charging process, the computation of the monomer den- 
sity by solving the modified diffusion equation would lead to the violation of the constraint 
given by Eq. ( IB-12P and make it difficult to deduce the energetic and entropic contributions 
to the free energy. Also, Eq. ( IB-12P is part of the reason that the dielectric constant of 
the medium has to be assumed to be independent of temperature, while carrying out this 
analysis. 

Having shown that the work done during the isothermal charging process can be used 
to compute the free energy of the system, it is clear that for the dual charging process 
considered here, they are related by F = Fg + -F^/2. Using this relation between the work 
done as a result of the charging and the free energy, the entropy can be calculated using Eq. 
( [B^ so that 
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Now, using the integral form for the Poisson equation obtained at the saddle point ip'{r) = 
—At[Ib f c?r'pg(r')/|r'|, the assumption that e is independent of T and w'pp = a/T, the above 
equations take the form 
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expressions for the entropic contributions can be cast in the form 

d (p'eir) 
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Carrying out integration by parts over A 
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Here, we have used the fact that the fields and charge density are zero for A = i.e., 
V^'(r) = Pe(r) = Wp(r) = for A = 0. Reahzing that the partition functon for a single small 
ion can be used to rewrite the entropy by 
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Similarly, the entropic contributions from the polymer can be rewritten using 
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Using Eqs. ( IB^ and ( IB^ 
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Carrying out the integrals over A, we get 
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where we have used the fact that Q'j = Q for A = 0. An important point to note here is 
that S turns out to be the entropy of the system relative to the ideal system (i.e., a system 
in the absence of interactions). In other words, F = E — TS where S = Stotai — Sid, as was 
already pointed out by Marcus in Ref.— . 

Now, it is easy to identify the entropic contributions arising from different components. 
Separating S into contributions from small ions and the polymer chain by writing 
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where quantities with "ic?" in the superscripts are the entropic contributions in the ideal 
system. In writing the entropy of small ions in terms of densities, we have used Eq. ( IB-26P 
after putting A = 1. Also, note that summing up the energy and entropy, the total free energy 
of the system obtained using the charging method described here differs from the free energy 
obtained within the saddle point approximation of SCFT by the entropic contributions of 
the ideal system in the absence of interactions, which is the most obvious choice for the 
reference frame during the computations of free energy in the field theory. In particular, the 
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free energy of the reference state comes out to be 
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APPENDIX C : Energy and entropy of single polyelectrolyte chain - incompressible 

system 

In Appendix B, we presented the derivation for the energetic and entropic contributions 
to the single polyelectrolyte chain system in the absence of solvent. Here, we present the 
generarlization of the technique to polyelectrolyte chain in the presence of solvent treated 
within the incomprcssibility constraint and point-like ions (i.e., Pp(r) + Ps(r) = po)- In this 
case, there are three excluded volume parameters Wpp, Wgs and Wps characterizing monomer- 
monomer, solvent-solvent and monomer-solvent interactions in contrast to just one in Ap- 
pendix B. Expressions for the electrostatic contributions remain the same and here, we focus 
on the contributions arising from the excluded volume interactions. 

Consider the two isothermal charging processes imagined in Appendix B i.e., the tra- 
ditional electric and excluded volume "charging", where charge of the ionic species and 
the excluded volume parameters characterizing the short range excluded volume effects are 
developed linearly and quadratically from to their final values in term of the charging 
parameter A. At any instance the excluded volume parameter between species i and j is 
w'j^j — X^Wij, where Wij is the final value of the excluded volume parameter. Also, the 
incomprcssibility constraint is forced at all instances during the charging process. 

Excluded volume energy for the system can be written as 



where Xpsb^ — — {wpp + Wss)/'2 is the Flory's chi parameter and it is taken to be of the 
form a/T from the assumed dependence of the excluded volume parameters Wij. 

Following the recipe presented in Appendix B, the effect of excluded volume interactions 
on the entropic contributions can be written as 
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where we have used w'^ir) = x'ps^^P'si^) + v' i'^) j '^'si^) = x'ps^^Ppi'^) + '7'(i") after adding 
^/^^^ ^{ppW+PsW} _ g ^]-^g jg^g^ step. Using the procedure described in Appendix B 
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where we have defined Q'^ = J rfr exp [— w^(r)] and p's{r) = rigexp [— Wg(r)]/Qg. An extra 
contribution due to the solvent appears in the expression for entropy so that 
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and the expressions for 
( IB-371) . respectively. 



Sl^^g and Spoiy — S^^^i are the same as in Eqs. ( IB-36|) and 
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FIGURE CAPTION 



Fig. 1.: Cartoons of the chain in the "free ends" (la) and "end- fixed" (lb) states. Spherical 
coordinates used to solve the SCFT equations are shown in Fig. Ic. 

Fig. 2.: (Color) Two dimensional monomer and electrostatic potential distribution for the 
single flexible polyelectrolyte chain in "free ends" [(a) and (6) , respectively ] and "end- 
fixed" state [(c) and (rf), respectively]. /^/^ = 3, a = 0.1, = O.IM, = 50, R/h = 4 
and Xps = 0.45. For plots (c) and (d), one end is anchored at [x, y] = [{R — 0.625)6, 0], 
which is shown by an arrow. 

Fig. 3.: (Color) Counterion and colon distribution for the single flexible polyelectrolyte 
chain in "free ends" [(a) and (6), respectively ] and "end-flxed" state [(c) and (d), 
respectively]. All the parameters are the same as in Fig. [2]i.e., Is/b = 3,a = 0.1, Cs = 
O.IM, = 50, R/b = 4 and Xps = 0.45. For the "end-flxed" state in plots (c) and (d), 
one end is anchored at [x, y] = [{R — 0.625)6, 0]. 

Fig. 4.: (Color) Monomer and electrostatic potential distribution for the single flexible 
polyelectrolyte chain in "free ends" [(a) and (6), respectively ] and "end-flxed" state 
[(c) and (d), respectively] at a higher monomer volume fraction compared to Fig. 
2. Parameters used to generate these plots are Is/b = 3,a = 0.1, c^ = O.IM, = 
100, R/b = 4 and Xps = 0.45. For plots (c) and (d), one end is anchored at [x,y] = 
[{R — 0.625)6, 0] (shown by arrow). 

Fig. 5.: (Color) Effect of conflnement on the monomer density distribution along x and y 
axes for the single flexible polyelectrolyte chain in "free ends" [(a) and (6) , respectively 
] and "end-flxed" state [(c) and (rf), respectively]. Parameters used to generate these 
plots are Is/b = 3,a = 0.1, c^ = O.IM, R/b = 4 and Xps = 0.45. For plots (c) and (d), 
one end is anchored at [x, y] = [{R — 0.625)6, 0]. 

Fig. 6.: Comparison of the monomer density distribution for the single flexible polyelec- 
trolyte chain in "free ends" and "end-flxed" states. Figs, (a) and (6) correspond 
to the density distribution along x and y axes, respectively for = 50. Similarly, 
Figs, (c) and (d) correspond to the density distribution along x and y axes, re- 
spectively for A^ = 100. All the other parameters are the same as in Fig. [2] i.e., 



Is/h = 3,a = 0.1, = 0.1M,i?/6 = 4 and Xps = 0.45. For the "end-fixed" state in 
plots, one end is anchored at [x, y] = [{R — 0.625)6, 0]. 



Fig. 7.: Comparison of the electrostatic potential (a), and the counterion {pc{x) + p+(x)) 
and colon (p_(x)) density distributions (6) for the single flexible polyelectrolyte chain 
in "free ends" and "end-fixed" state, respectively. All the parameters are the same as 
in Fig. li.e., h/b = 3, a = 0.1, c, = 0.1M,R/b = 4, iV = 100 and Xps = 0.45. 

Fig. 8.: (Color) Effect of and R on free energy barriers for the chain end to be localized on 
the surface of a neutral spherical cavity. Is/b = 3, a = 0.1, Cg = O.liVf and Xps = 

0.45. 

Fig. 9.: (Color) Dominance of conformational entropy to free energy barriers. For these 
figures, Is/b = 3, a = 0.1, = O.IM, = 0.45, and Figs, (a) and (b) correspond to 
R/b = 4 and R/b = 6, respectively. The net free energy barriers in the figures are the 
same as in Fig. [HI 

Fig. 10.: (Color) Dependence of the free energy barriers on the degree of ionization of 
the polyelectrolyte chain. Parameters used to obtain these plots are: Is/b = 3, = 
O.IM, R/b = 4:,Xps = 0A5. 

Fig. 11.: (Color) Effect of the added salt concentration on the free energy barriers. Pa- 
rameters used to obtain these plots are: Is/b = 3,a = 0.1, R/b = 4, Xps = 0.45. 

Fig. 12.: (Color) Dependence of the free energy barriers on Bjerrum length, which char- 
acterizes the electrostatic interaction strength between charged species. Parameters 
used to obtain these plots are: a = 0.3, Cg = O.IM, R/b = 4, Xps = 0.45. 

Fig. 13.: (Color) Effect of Flory's chi parameter (characterizing the exlcuded volume inter- 
actions between monomers and solvent molecules) on the free energy barriers is shown 
here. Parameters used to obtain these plots are: Is/b = 3,a = 0.1, Cg = O.IM, R/b = 
4. 




FIG. 1: Cartoons of the chain in the "free ends" (la) and "end-fixed" (lb) states. Spherical 
coordinates used to solve the SOFT equations are shown in Fig. Ic. 




FIG. 2: (Color) Two dimensional monomer and electrostatic potential distribution for the single 
flexible polyelectrolyte chain in "free ends" [(a) and (b), respectively ] and "end- fixed" state [(c) 
and (d), respectively]. Is/b = 3,a = 0.1, c<j = O.IM, iV = 50,R/b = 4 and Xps = 0.45. For plots 
(c) and (d), one end is anchored at [x,y] = [{R — 0.625)6,0], which is shown by an arrow. 
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FIG. 3: (Color) Counterion and coion distribution for the single flexible polyelectrolyte chain in 
"free ends" [(a) and (6), respectively ] and "end-fixed" state [(c) and (d), respectively]. All the 
parameters are the same as in Fig. [2]i.e., Is/b = 3, a = 0.1, Cg = O.lAf, N = 50, R/b = 4 and Xps = 
0.45. For the "end-fixed" state in plots (c) and (d), one end is anchored at [x, y] = [{R — 0.625)6, 0] . 




FIG. 4: (Color) Monomer and electrostatic potential distribution for the single flexible poly- 
electrolyte chain in "free ends" [(a) and (6), respectively ] and "end-fixed" state [(c) and (d), 
respectively] at a higher monomer volume fraction compared to Fig. 2. Parameters used to gen- 
erate these plots are Is/b = 3,a = 0.1, = O.IM, N = 100, R/b = 4 and Xps = 0.45. For plots (c) 
and (d), one end is anchored at [x,y] = [{R — 0.625)6,0] (shown by arrow). 
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FIG. 5: (Color) Effect of confinement on tlie monomer density distribution along x and y axes 
for the single flexible polyelectrolyte chain in "free ends" [(a) and (b), respectively ] and "end- 
fixed" state [(c) and (d), respectively]. Parameters used to generate these plots are Is/b = 3, a = 
0.1, Cs = O.IM, R/b = 4 and Xps = 0.45. For plots (c) and (d), one end is anchored at [x,y] = 
[(i?- 0.625)6,0]. 
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FIG. 6: Comparison of the monomer density distribution for the single flexible polyelectrolyte 
chain in "free ends" and "end-fixed" states. Figs, (a) and (6) correspond to the density distribution 
along X and y axes, respectively for = 50. Similarly, Figs, (c) and (d) correspond to the density 
distribution along x and y axes, respectively for N = 100. All the other parameters are the same 
as in Fig. [2]i.e., Is/b = 3,a = 0.1, = 0.lM,R/b = 4 and Xps = 0.45. For the "end-fixed" state 
in plots, one end is anchored at [x, y] = [{R — 0.625)&, 0]. 
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FIG. 7: Comparison of the electrostatic potential (a), and the counterion {pc{x) + and 
colon density distributions (6) for the single flexible polyelectrolyte chain in "free ends" 

and "end-flxed" state, respectively. All the parameters are the same as in Fig. |3]i-e., Is/h = 
3, a = 0.1, Cs = O.IM, R/h = 4, = 100 and Xps = 0.45. 
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FIG. 8: (Color) Effect of N and R on free energy barriers for the chain end to be localized on the 
surface of a neutral spherical cavity. Is/b = 3, a = 0.1, Cg = O.IM and Xps = 0.45. 
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FIG. 9: (Color) Dominance of conformational entropy to free energy barriers. For these figures, 
Is/b = 3, a = 0.1, c<j = 0.1M,Xps = 0.45, and Figs, (a) and (b) correspond to R/b = 4 and 
R/b = 6, respectively. The net free energy barriers in the figures are the same as in Fig. [51 




FIG. 11: (Color) Effect of the added salt concentration on the free energy barriers 
used to obtain these plots are: Is/b = 3,a = 0.1, R/b = 4, Xps = 0.45. 
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FIG. 12: (Color) Dependence of the free energy barriers on Bjerrum length, which characterizes 
the electrostatic interaction strength between charged species. Parameters used to obtain these 
plots are: a = 0.3, = 0.1M,i?/6 = A.Xps = 0.45. 
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FIG. 13: (Color) Effect of Flory's chi parameter (characterizing the exlcuded volume interactions 
between monomers and solvent molecules) on the free energy barriers is shown here. Parameters 
used to obtain these plots are: Is/b = 3, a = 0.1, Cg = 0.1M,R/b = 4. 



